Macroscopic current fluctuations in stochastic lattice gases 
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We study current fluctuations in lattice gases in the macroscopic limit extending the dynamic 
approach to density fluctuations developed in previous articles. More precisely, we derive large de- 
viation estimates for the space-time fluctuations of the empirical current which include the previous 
results. Large time asymptotic estimates for the fluctuations of the time average of the current, 
recently established by Bodineau and Derrida, can be derived in a more general setting. There are 
models where we have to modify their estimates and some explicit examples are introduced. 
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In a series of recent papers 0, S 01 we have devel- 
oped a dynamical approach to macroscopic fluctuations 
of thermodynamic variables in non equiUbrium steady 
states of stochastic lattice gases. The main object stud- 
ied in these papers is the empirical density for which we 
established a large deviation principle. From this prin- 
ciple we obtained several results: we have shown that 
the entropy functional satisfies a Hamilton-Jacobi equa- 
tion, we have extended Onsager-Machlup theory and 
the minimum dissipation principle j^] to stationary non 
equiUbrium states. 

Another observable of great physical interest is the 
macroscopic current flowing through the system |^ 0, 
H, In the present article we develop, in the same 

spirit of 0, Q , a large deviation principle for the empiri- 
cal current and discuss some of its consequences. We dis- 
cuss the asymptotic probability, as the number of degrees 
of freedom increases, of observing a space-time depen- 
dent current profile different from the typical value. This 
new principle includes the previous one for the density 
and leads to a unified approach to fiuctuations in steady 
states both equilibrium and non equilibrium. From this 
asymptotics we derive a more general form of a large 
fiuctuati on p rinciple recently obtained by Bodineau and 
Derrida [llj for the large time fluctuations of the time 
averaged current. We show that there are models where 
one has to replace the large deviation functional obtained 
in ^3 1 that we denote by U{J), by its convex envelope 
U**{J). We discuss some examples where U** < U; this 
is naturally interpreted as a dynamical phase transition. 
However more complex situations are also conceivable. 

In a forthcoming more detailed paper we shall 
study also the behavior of the current large deviation 
functional under time reversal and connect it to the 
well known fluctuation theorem for entropy production 
of Gallavotti and Cohen El [1113. 

As the basic microscopic model we consider a stochas- 
tic lattice gas with a weak external field and particle 



reservoirs at the boundary. More precisely, let A C M'* 
be a smooth domain and set A^r = iYAnZ'*; we consider 
a Markov process on the state space X^" , where A is a 
subset of N, e.g. X = {0, 1} when an exclusion principle 
is imposed. The number of particles at the site x £ A^v 
is denoted hy r]x € X and the whole configuration by 
7] G A'^" . The dynamics of the particles is described by a 
continuous time Markov process on the state space A^™ 
with transition rates Cx,y{ri) from a configuration ry to 
a'^'yrj that is the configuration obtained from rj by mov- 
ing a particle from x to a nearest neighbor site y. Similar 
rates describe the appearance or loss of a particle at 
the boundary site x. The reservoirs are characterized by 
a chemical potential 7. We assume that the rates sat- 
isfy the local detailed balance condition Ts^ with respect 
to a Gibbs measure associated to some Hamiltonian 7i. 
Typically, for a non equilibrium model, we can consider 
A the cube of side one and the system under a constant 
force E/N. Moreover we choose the chemical potential 
7 so that jiy/N) = 70 if the first coordinate of y is 0, 
j{y/N) = 71 if the first coordinate of y is N, and impose 
periodic boundary conditions in the other directions of 
A. 

The fluctuation theory of stochastic lattice g ases, as 
discussed below, is expected to apply to a wider class 
of non equilibrium systems with conservation laws. For 
instance, the Hamiltonian system discussed in 16] is re- 
duced to the analysis of a Markov process. We introduce 
the empirical density tt^ associated to a microscopic con- 
flguration ?] G A'^" by requiring for each smooth function 
G : A ^ M, 

{n^,G) = I dun^{u)Giu) = ^ ^ Gi^/N)n. 



xeAf, 



so that TT^ {u) is the local density at the macroscopic 
point u = x/N in A. 

Consider a sequence of initial configurations rj^ such 
that 7r^(?7^) converges to some density profile po- Under 
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diffusive scaling the empirical density at time t converges, 
as N —^ oo, to p — p{t, u) which is the solution of 



dtp = y 



x{p)VV ^V{p) (1) 



with initial condition pQ. Here D is the diffusion ma- 
trix, given by the Green-Kubo formula, see [H II.2.2], 
X is the conductivity, obtained by linear response the- 
ory, see 0, II. 2. 5], and W the external field. We em- 
phasize that these transport coefficient are defined in 
terms of the equilibrium Gibbs measure. In particular 
if we denote by So{p) the entropy associated to H, the 
usual Einstein relation D{p) = R~^{p)x{p) holds, here 
R{p) = Sq{p)~^ is the compressibility. The interaction 
with the reservoirs appears as a boundary condition to 
be imposed on solutions of (0). More precisely, we re- 
quire that S'q{p{u)) = 7(m), u £ OA, here i9A denotes the 
boundary of A and we recall that 7 is the chemical po- 
tential of the reservoirs. The non equilibrium stationary 
profile p is the unique stationary solution of 

Denote by p^ the non equilibrium stationary ensem- 
ble, i.e. the invariant measure of the microscopic process. 
The probability of density fluctuations can be described, 
for large N, by an entropy functional S{p) which depends 
on the local density of particles p. The stationary pro- 
flle p corresponds to a critical point of S so that it is a 
minimum of the entropy; here we use the probabilist sign 
convention for S which is opposite to physicist conven- 
tion. Hereafter, we normalize the entropy S{p) so that 
Si/5) = 0. As shown in 0,0, the entropy functional S 
satisfies the following Hamilton-Jacobi functional deriva- 
tive equation 



1 /^SS , ,^5S\ 







(2) 



where (•, •) means integration with respect to the space 
variable u £ A. It shows that once we know the hydro- 
dynamic equation, we can obtain a closed macroscopic 
description. 

The probability to observe a macroscopic trajectory 
different from the hydrodynamic behavior ^ is expo- 
nentially small in N"^ and given by 



p, [0,T]) ^exp{ -iV%,T](p)} (3) 



where « denotes closeness in some metric, ~ logarithmic 
equivalence as TV 00, and stands for the distri- 
bution of the stationary process, namely we consider an 
initial condition distributed according to the invariant 
measure p^ . The rate functional /[o.t](p) is given by 



/[o,T](p) = 5(p(0)) + - dt{VH,x{pWH) (4) 

where the external potential H has to chosen so that p 
solves 



dtp = V{p) - V • (x(p)Vi7) 



with the same boundary conditions as the hydrodynamic 
equation Equations lEJ^© represent a dynami- 

cal generalization of Einstein formula for thermodynamic 
fluctuations. 

We give now a parallel discussion for the current. De- 
note by A/^'^ the number of particles that jumped from x 
to y in the macroscopic time interval [0, t]. Here we adopt 
the convention that M^''^ represents the number of par- 
ticles created at y due to the reservoir at x if a; ^ Aat, 
y G An and that A/'/''*' represents the number of particles 
that left the system at x by jumping to y if a:; G Ajv, 
y ^ An- The difference J^'^ — N^'^ — M^'^ represents 
the total current across the bond {x, y} in the time in- 
terval [0, t]. In other words, given a path 77(5), < s < t, 
the instantaneous current dJ^ '^ /dt is a sum of (5-function 
localized at the jump times across the bond {x,y} with 
weight respectively —1, if a particle jumped from x 
to y, respectively from y to x. 

Fix a macroscopic time T and denote by J''^ the em- 
pirical measure on [0, T] x A associated to the current. 
For smooth vector fields G = (Gi, . . . , Gd), the integral 
of G with respect to J'^ , denoted by {G)^ is given by 

J^{G) ^ [ dt [ duG{t,u) ■ J'^{t,u) 

Jo J A 

cT 



See / dtGs,x/N)'''' 
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dt 



where • is the inner product in R'', is the canonical 
basis, and we sum over all x such that either x G A^v or 
x + 6; G An- We normalized so that it is finite as 

N ^ OQ. 

Given a density profile p let us denote by 



J{p)^-\d{p)^p + x{pWv 



(6) 



the current associated to p- The hydrodynamic equation 
can then be written ds, dtp + V ■ J [p] — Q . If we con- 
sider an initial conflguration rj^ such that the empirical 
density ir^ [rj^) converges to some density profile po and 
denote by p{t) the solution of then the empirical cur- 
rent {t) converges, as iV ^ cxo, to J{p{t)), the current 
associated to the solution of the hydrodynamic equation 
([T|). If we let the macroscopic time diverge, t —^ 00, we 
have J (pit)) J{p). 

We next discuss the large deviation properties of the 
empirical current. Fix a smooth vector field j : [0,T] x 
A ^ M'' and a sequence of configurations 77^ whose em- 
pirical density converges to some profile po- The large 
deviation principle for the current states that 

P^4j''it,u)^jit,u)) ^expi-N'j^o.rM (7) 
where the rate function is given by 



(5) ^[0,T]ij) 



dt{[j^J{p)],xipr'[j-Jip)]) (8) 
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in which J{p) is given by © and p — p(t, u) is obtained 
by solving the continuity equation dtp + V ■ j — with 
initial condition p(0) = Pq. Of course there are com- 
patibility conditions to be satisfied, for instance if we 
have chosen a j such that p{t) becomes negative for some 
te[0,r] thenZ[o,T](i) = +00. 

We present here a heuristic derivation of 131-© • Fix 
a current j; in order to make j typical, we introduce an 
external field F. Let p be the solution of 



dtp + \7-j = 
p{Q,u) = pa{u) 

and i^" : [0, T] X A ^ R'' be the vector field such that 

J = J{p)+x{p)F 



(9) 



= --D{p)Vp + x{p){'^V + F} 



We introduce a perturbed measure 



which is ob- 



tained by modifying the rates as follows 

Following a similar argument as the one for the large 
deviation principle of the empirical density 0, one can 
show that 



aN 



.N,F 



exp 



{ 



dt 



{F,x{p)F)] 



= exp{-7V%.T](j)} 



Moreover, under 
Therefore, 



T>N,F 



as iV — > 00, converges to j. 



We emphasize that now we allow non-gradient external 
fields F, while in the large deviation principle for the 
empirical densi ty it is sufficient to consider gradient 
external fields |l8l ll9| . The latter is therefore a special 
case and can be recovered from (|7|)-(|Hll- 

We want to study the fluctuations of the time average 
of the empirical current over a large time interval [0, T]; 
the corresponding probability can be obtained from the 
space-time large deviation principle Q. Fix some di- 
vergence free vector field J = J{u) constant in time 
and denote by At,.j the set of all currents j such that 
dt j{t,u) = J{u). The condition of vanishing di- 
vergence on J is required by the local conservation of the 
number of particles. By the large deviations principle 
{Tj), for T large we have 



P^«(^/ dtj^{t) 



J) -exp{ -iV'*T$(J)} (10) 



where the logarithmic equivalence is understood by send- 
ing first N ^ 00 and then T — > 00 . In |l2l | we show that 
for the so-called zero range process the limits can be 
taken in the opposite order; we expect this to be true in 
general. The functional $ is given by 



<i>(J) = lim inf 



(11) 



By a standard sub-additivity argument it is indeed easy 
to show the limit exists. We now prove that $ is a con- 
vex functional. Let J = pJi -I- (1 — p)J2- we want to 
show that $(J) < p$(Ji) + (1 - p)$(J2). Let us caU 
t G [0,pT] (respectively (j2(t), P2(t)), t € 
[0,(1 — p)T]) the optimal path of current and density 
which implements the inf in (|ll|l at Ji (resp. J2). We 
then consider a path (j(t), p{t)) which spends a time in- 
terval pT following (ji,pi) and a time interval (1 — p)T 
following (^2,^2) (and a finite time to go continuously 
from (ji, pi) to (j2, P2))- With such a path we get 



1 



+ y2.[o.(i_j, 



)T] (j2) 



< p$(Ji) + (l-J5)$(J2 



In the last step we used the existence of the limit Hll|l . 

We next study the variational problem on the right 
hand side of l|ll|l . We begin by deriving an upper bound. 
Given p = p{u) and J = J(m), V • J = 0, let us introduce 
the functionals 

u{p.J) = \{J - J{p)Mp)-\J - J{p)]) (12) 



U{J) = inf U{p, J) 
p 



(13) 



where the minimum in (|13|l is carried over all profiles p 
satisfying the boundary conditions and J{p) is given by 
©. When J is constant, that is in the one-dimensional 
case, the functional U is the one introduced in 11]. 

By choosing a suitable path j{t,u) G At^j we first 
show that 



$( J) < C/(J) 



(14) 



The strategy is quite simple. Let p — p{J) be the den- 
sity profile which minimizes the variational problem l|13|l . 
Given the initial density profile po, we construct a path 
j = j{t, u), {t, u) G [0, r] X A as foUows 



J 



Jit) 



T-2t 



if < i < T 
J if T <t <T -T 

j a T -T<t<T 



where j is a vector field such that rV- j — po—p and r > 
is some fixed time. It is now straightforward to verify 
that j G At,j, as well as limT-^oo ^ ^[o,T]{j) = U{J). 

From H14|l and the convexity of ^(J) it immediately 
follows that 



$(J) < U**{J) 



(15) 
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We next discuss a lower bound for the variational prob- 
lem Hll() . We denote by U and U the same functionals as 
in H12|l - (|13() . but now defined on the space of all currents 
without the conditions of vanishing divergence. Let also 
U** be the convex envelope of U. 

Let j e At. J, by the convexity of U** in the set of all 
currents, we get 



1 



2^[0,T](j) 



1 

T 



dtU{p{t),j{t)) > 



dt U{j{t)) 



(16) 



The upper and lower bounds H15() and (|16|l are, in 
general, different. For a divergence free J we have 
U{J) = U{J) but since the convex envelopes are consid- 
ered in different spaces, we only have U**{J) < U**[J). 

To understand the physical meaning of the convex 
envelope in l(T^ . suppose J ~ pJi + (1 — p)J2 and 
U{J) > U**{J) = pU{.h) -f (1 - p)f/(J2) for some 
P, J\tJi- The values p, Ji, J2 are determined by J and U . 
In addition we assume that U** = [/**. If we condition 
on observing an average current J, the corresponding 
density profile is not determined, but rather we observe 
with probability p the profile p{Ji) and with probability 
1 — p the profile p{J2)- When U is not convex we have 
thus a situation in which the time averaged current J is 
realized with the coexistence of two dynamical regimes: 
we have a dynamical phase transition. 

The derivation of the upper bound shows that, if U is 
not convex, our result differs from the one in On 
the other hand if U**{J) < U**{J) it is possible that one 
can improve it by exploring currents with non vanishing 
divergence. In such a situation it is not clear to us if <& 
can be directly related to U. 

We can consider the large time behavior of the em- 
pirical current as in equation H1U|) with the additional 
constraint that the associated density is asymptotically 
time independent. In such a case it is not difficult to show 
that (|10|l holds with ^ = U. With this extra constraint 
we are in fact forbidding the system from oscillating. 

In the models where the diffusion coefficient D{p) is 
constant and the mobility xip) is concave, for example 
in the symmetric simple exclusion where x = p(l — p), 
it is not difficult to see that U is convex. Therefore in 
these cases ^ — U. In we shall show that in the 
Kipnis-Marchioro-Presutti model [i^, 0| at equilibrium 
U{J) is convex while U{J) is not. 

We next discuss an example, with a non-concave Xj 
where the functional U is not convex. We fix equilibrium 
boundary conditions po = Pi = P- We take D{p) = 1 
and x{p) ^ smooth function with x(0) — = which 



has a non concave part x{p) = Ke~''^P, for p in a given 
interval, where C is a positive parameter. An explicit 
calculation gives that 



U{J) = ^e^^FiCJ) 



(17) 



where F{z) — z^^^, if J is in an appropriate interval. 
Since the second derivative of F{z) can be negative, U{J) 
is not convex. 
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